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ABSTRACT 

We use galaxy catalogues constructed by combining high-resolution N-body simulations with 
semi-analytic models of galaxy formation to study the properties of Long Gamma-Ray Burst 
(LGRB) host galaxies. We assume that LGRBs originate from the death of massive young 
stars and analyse how results are affected by different metallicity constraints on the progenitor 
stars. As expected, the host sample with no metallicity restriction on the progenitor stars 
provides a perfect tracer of the cosmic star formation history. When LGRBs are required 
to be generated by low-metallicity stars, they trace a decreasing fraction of the cosmic star 
formation rate at lower redshift, as a consequence of the global increase in metallicity. We 
study the properties of host galaxies up to high redshift (~ 9), finding that they typically have 
low-metallicity (Z < O.5Z ) and that they are small (M < 1O 9 M ), bluer and younger 
than the average galaxy population, in agreement with observational data. They are also less 
clustered than typical galaxies in the Universe, and their descendents are massive, red and 
reside in groups of galaxies with halo mass between 10 13 M Q to 10 14 M Q . 

Key words: gamma-rays: bursts - host galaxies . 



1 INTRODUCTION 

Gamma-ray b ursts (GRBs) are the mo st energetic explosions in 
the Universe (Zhang & Meszaros 2004). As such, they offer ex- 
citing possibilities to study astrophysics in extreme conditions, 
e.g., r adiative processes in highly relativistic ejecta dHuang et al.l 
l200d : lFan & Piranll2008l . and references therein). Because of their 
very large luminosity, GRBs represent "cos mological" events , 
which have been det ected up to z ~ 8.2 l lTanvir et alj 1 20091 ; 
ISalvaterra et al.l 12009). It has been proposed that some tight cor- 
relations among GRB parameters can make them "standard can- 
dles" for probing the Universe to the high-redshift regime that 
supernovae la can not attain (e.g. Ghirlanda et al. 20 041; [Dai et alj 
20041 Zhand2007l and references therein. See however lLil2007jbl 



2008a) 



It is well-known that the d uration distribution of GRBs is bi- 
modal jKouveliotou et alJI 19931) . dividing GRBs into two classes: 
long GRBs (hereafter LGRBs) and short GRBs, depending on 
whether their durations are longer or shorter than a few sec- 
onds. The observed properties of host galaxies of short and long 
GRBs indicate that they have different progenitors. LGRBs are 
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typically found in st ar-forming galaxies, predominantly i rregu- 
lar dwarf galaxies (jConselice et all 120051 : iFruchter et alj 120061 : 
IWainwright et alj|2007l) . In contrast, short GRBs are found in both 
early-type and late-type galaxies. Many models have been pro- 
posed for explaining the origin of these two classes of GRBs. The 
currently favourite hypotheses are that short GRBs are produced 
by the merger of compact objects - betwee n two neutron stars or 
between a neutron star and a black hole jLi&Paczvrtskilll998l : 
lO'Shaughnessv et al.l2 008). while LGRBs originate from the death 
of massive stars (with low metallicity), such as Wolf-Rayet stars. 
This scenario for the fo rmation of LGRBs is usually referred to as 
the "c ollapsar model" jYoon et alj|2006l 120081 ; IWooslev & Hegerl 
120061) . 

Observational data are consistent with the hypothesis of a 
LGRB -supernova connection: at le ast some LGRBs a re associ- 
ated with core-co llapse supernovae jGalamaetal.ll 19981 ; 10 120061 : 
IWooslev"&^ Heger 2006, and references therein). In addition, all su- 
pernovae associated with GRBs are Type Ic, which supports the hy- 
pothesis of Wolf-Rayet stars as progenitors of LGRBs. Because of 
their connection with supernovae, LGRBs are potential tracers of 
the cosmic star formation history (Totani|[l997l ; IWiiers et alj[l998l : 
iMao & Molll998l ; |Porciani & Madaull200ll) . To date, there are 130 
GRBs w ith known redshift an d 50 with estimated host galaxy stel- 
lar mass JSavaglio et ai]2 008). Given the difficulty in detecting and 
localizing short GRBs, most of the observational studies about host 
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WMAPl 


WMAP3 




0.25 


0.226 


n A 


0.75 


0.774 


n b 


0.045 


0.04 


fg 


0.9 


0.722 


ft 


0.73 


0.743 


n 


1 


0.947 



Table 1. Cosmo logical parameters of the two simulations used in 
Wang et al. (2008). Q m , Ha, represent the density of matter, dark en- 
ergy, and baryons respectively, erg and n are the amplitude of the mass den- 
sity fluctuations, and the slope of the initial power spectrum. The Hubble 
constant is parameterised as Hq = 100 h kms _1 Mpc _1 . 



galaxies are for LGRBs, which will be the focus of this work. Stud- 
ies of the physical properties of GRBs are not easy as they require 
deep targeted observations at high redshift. In addition, the proba- 
bility for a chance superposition of GRBs and galaxies on the sky 
is significant fo r high-z GRBs, and ~ 3% for galaxies at z < 1.5 
dCampisi & OI2008I) . 

The observational information gathered so far indicates that 
most LGRBs are found in faint star forming galaxies domi- 
nated by young stellar populations with a sub-solar gas-phase 
metallicities, although there are a few host galaxies with higher 
metal content ( Prochaska et al. 2004; Wolf & Podsiadlowski 2007; 
Fvnbo et al.ll2006l:|Price et a l. 2007; Savaglio et al. 2003; S avaglid 



2006; Savaglio et al. 2008; Stanek 



Savaglio et a. 
: et all l2006t 



and references 
therein). 

In this work, we analyse the properties of host galaxies 
of LGRBs, using a galaxy catalogue constructed by combining 
high-resolution N-body simulations with a semi-analytic model of 
galaxy formation. In particular, we use the models discussed in 
Wang et al. (2008) for two cosmological models with parameters 
taken from the first-year an d the third-year Wil kinson Microwave 
Anisotropy Probe (WMAP) ( Spergel e t alj2003h measurements. To 
select candidate host galaxies of LGRBs, we extract from the avail- 
able semi-analytic galaxy catalogues the information for the age 
and metallicity of newly formed stars and we adopt the collap- 
sar model. We built three samples of host galaxies with different 
metallicity thresholds, and we compare the properties of the se- 
lected galaxie s with observational data, in particular the data in 
I Savaglio et alj |2008). Compared with previous theoretical studies, 
the simulations used in our study have the largest volume, and they 
also allow us to explore the cosmological dependence. Finally, the 
information available from the semi-analytic catalogues enable us 
to study the clustering and descendent properties of LGRB hosts. 

The paper is organised as follows. We present in section 2 
the simulated galaxy catalogues used in this work. In section 3, 
we describe the method for selecting long GRB host galaxies. We 
describe our results in section 4. We discuss our results and give 
our conclusions in section 5. 



2 THE SIMULATED GALAXY CATALOGUES 

In this study, w e use the galaxy catalogues constructed by 
IWang et alj ( 2008) for two simulations with cosmological param- 
eters from the first and third-year WMAP results. The two sets 
of cosmologic a l para meters are listed in Table [T] As discussed 
in IWang et all |2008), the most significant differences between 
WMAPl and WMAP3 cosmological parameters are a lower value 



of <rg and a redder (smaller) primordial power spectrum index n 
in WMAP3, resulting in a significant delay for structure formation. 
Both simulations correspond to a box of 125 ft _1 Mpc comoving 
length and a particle mass 8.6 x 10 s (WMAPl) and 7.8 x 10 s M 
(WMAP3). The softening length is 5 h^ 1 kpc in both simulations 
(see Table 2 in Wang et al. 2008). Simulation data were stored 
in 64 outputs, that are approximately logarithmically spaced in 
time between z — 20 and z = 1, and linearly spaced in time 
for z < 1. Each simulation output was analysed with the post- 
process ing software original ly developed for the Millennium Sim- 
ulation dSprrngel et alj|2005h . 

Merging history trees for self-bound structures extracted from 
the simulations were used as input for the Munich semi-analytic 
model of ga l axy f o rmation. Interested r eaders are referred to 
ICroton et all d2006t) . IDe Lucia & Blaizotl d2007l) and references 
therein for details on the physical processes explicitly modelled. 
Previous work has shown that the galaxy population predicted 
by this particular model provides a reasonably good match with 
the observed local galaxies prope rties and relations among stel 



lar m ass, gas mass, and metallicity jPe Lucia. Kau ffmann & White 



2004), luminosity, co lour, morphology distributions ( ICroton et al 



2006; De Lucia et aT]|2006l). and the observed two-p oint correla- 
tion functions dSpringel et al]2005l ; IWang et alj2 008). In addition, 
Kitzbichle r~& White! d2007l) have shown that the model also agrees 
reasonably well with the observed galaxy luminosity and mass 

function at higher redshift. 

W e remind the reader that the models discussed in Wang et al. 
(2008) adopt the same physics, but different combinations of model 
parameters are used for the simulations with WMAPl and WMAP3 
cosmology. The WMAPl simul ation uses the same param eters 
(and physical model) adopted in IDe Lucia & Blaizotl d2007l) . For 
the WMAP3 simulation, the model adopts lower supernovae and 
AGN feedback efficiencies in order to compensate for the delay 
in structure formation obtained with a lower erg. This combination 
of model paramet ers corresponds to the WMAP3B model used in 
IWang et alj fc008h . The alternative model WMAP3C used in that 
paper leads to very similar results and so will not be discussed fur- 
ther in this paper. In the following, we limit our analysis to galaxies 
with stellar mass larger than 2 x 10 8 M©, which is above the reso- 
lution limit of the N-body simulations used. 

We note that the most recent cosmological model from the 
five-year data of WMAP is between WMAPl and WMAP3, so the 
results from the two simulations used here are expected to bracket 
results from a simulation with the 5-year WMAP cosmology. 



3 IDENTIFICATION OF LGRB HOST GALAXIES 

In order to identify candidate GRB host galaxies, we adopt the 
collapsar model for LGRBs: all young stars with mass > 30 Mq 
ending their life with a supernova should be able to create a BH 
remnant. If the collapsar has high angular mom entum, the forma- 
tion o f the BH is accompanied by a GRB event dYoon et alfeood 
2008). As mentioned in Sec. Q] recent studies on the final evolu- 
tionary stages of massive stars have suggested that a Wolf-Rayet 
(WR) star can produce a LGRB if its mass loss rate is small, which 
is possible only if the metallicity of the star is very low. When 
metallicities are lower than ~ 0.1 — 0.3 Zq, the specific angular 
momentum of the progenitor allows the loss of the hydrogen en- 
velope while preserving the heliu m core dWooslev & Hegerll2006l: 
iFrver. Wooslev & Hartmanrl 1 1 9991) . The loss of the envelope re- 
duces the material that the jet needs to cross in order to escape, 
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JVIMg 1 ] m,, m in 




SNe 


7.421 X 10~ 3 8 


100 


BH 


1.035 x 10" 3 30 


100 



Table 2. Number of SNe and BHs per solar mass of stars formed, using a 
Salpeter IMF. 



while the helium core should be massive enough to collapse and 
power a GRB. 

In order to construct our host galaxy sample, we have ex- 
tracted from the available semi-analytic catalogues the information 
about the age and metallicity of all stars. We have then created the 
following host samples: 



(i) HOST1, obtained by selecting galaxies containing stars with 
age < t c = 5 x 10 7 yr; 

(ii) HOST2, including galaxies with stars of age < t c and metal- 
licity Z sC 0.3Z Q ; 

(iii) HOST3, defined by selecting galaxies containing stars with 
age < t c and metallicity Z ^ 0.1Z@. 

In order to count the number of GRB events in each galaxy, 
we make use of two important pieces of information: (a) the 
rate of GRB with respect to the SNe explosions (without any 
cut in metallicity for the progenitors stars); (b) the rate of very 
massive stars (producing remnant BHs) as a function of redshift 
and metallicity with respect to the total number of SNe events. In 
this way we are able to count how many BHs with low metallicity 
progenitors will produce GRBs. 

We assume a Salpeter Initial Mass Function (IMF) and com- 
pute the number of stars ending their lives as supernovae (SNe) or 
as black-holes (BHs) per unit mass of stars formed by: 



N 



(m») dm* 



(1) 



m*(b(m t ) dm, 



where 4>(m*) is the IMF, m* jm i n is the minimum initial mass to 
form a supernova or a black hole, and ?Tl* , max is the upper limit of 
the mass f unction. We take rru.min to be 8 Mq for SNe and 30 Mq 
for BHs dFrveretal.lll999l) . and m,, max = 100 Mq (the lower 
limit of the IMF is taken to be 0. IMq). Eq.Q] evaluated for progen- 
itor stars of all metallicities, provides the numbers listed in Table[2] 
Without any restriction on the metallicity of the progenitor stars, 
the relative number of BHs and supernovae is ~ 14 per cent (140 
BHs per 1000 SNe). When considering the meta llicity threshold for 
BHs, the rate is reduced and varies w ith redshift dLanger & Normanl 
l2006l:IWolf & Podsiadlowsirill2007l) . This is illustrated in Fig.Q]for 
the two different metallicity thresholds adopted in our study. At 
low redshift, there are from 0.1 to 6 low-metallicity BHs formed 
per 1000 SNe. The rate grows with increasing redshift, reaching 
values between ~ 10 (Z < 0.1 Zq) and ~ 70 (Z < 0.3 Zq) at 
redshift ~ 9. 

Integrating over the redshift, Fig. Q] allows us to calculate the 
average ratio between the number of black holes from different 
progenitors and the number of supernovae in the Universe. How- 
ever, not every black hole will produce a LGRB. To normalise the 
LGRB abundance, we assume that the rate of GRB per SNe is on 
average (over all cosmic times) of about 1 GRB event every 1000 
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Figure 1. Number of low-Z BHs per 1000 SNe as a function of redshift, for 
progenitor stars with metallicity lower than . 3 Zq and . 1 Zq. Results are 
shown for the model WMAP3B. Similar results are obtained in the WMAP1 



SNe dPorciani & Madauj200ll;lLanger & Normarj200f3) . The num- 
ber of LGRBs relative to the black holes in the three samples (with 
or without the metallicity cutoff) is given by 



Rgkb = 



GRBs 
BHsi ow z 



GRBs 
lOOOSNe 



x ^ lOOOSNe^ 
BHsiowz 



(2) 



where the last term is obtaining by integrating over redshift (or 
cosmic time). For the HOST1 sample, the last term is a constant 
(1/140), and thus Rgrb ~ 0.007; for the other two samples, we 
obtain 7?grb ~ 0.056 and ~ 1. 

For each galaxy in the simulation box, we can count how many 
BHs are produced from low metallicity progenitors and then obtain 
the corresponding number of LGRBs: 

A^lgrbs = BHsi ow z x -Rgrb- 

It is important to note that our model provides only an up- 
per limit to the number of LGRB events because the formation of 
BHs and low metallicity are only two of many requirements for 
the production of LGRBs. E.g., high spin of the progenitor star is 
another requirement that cannot be handled in our model. In the 
following, we consider each galaxy hosting at least one GRB event 
(A^lgrbs > 1) as a 'host galaxy'. 

There are relatively few galaxies with A^lgrbs < 1 and their in- 
clusion does not change our results significantly. 
We stress that we are not using the (average) galaxy metallicity to 
select our host galaxy samples, but the metallicity of each 'pocket' 
of stars formed at each time-step. Stars are assumed to form with 
the metallicity of the interstellar medium at the time of star for- 
mation, and the model adopts an instantaneous recycling approxi- 
mation for metal enrichment. So, the gas-phase metallicity of host 
galaxies will generally be higher than the metallicity threshold we 
have adopted for our samples. 

Note that the LGRB rate computed above is not directly compara- 
ble to the observed rate because that would require us to take into 
account many unknown factors lik e the jet angle, and to i nclude 
any possible observational bias (see lLapi et al. 2008; Li 2008bT). 
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4 RESULTS 

In this section, we discuss the physical properties of LGRB hosts 
selected using the procedure described in Sec. [3] In particular, we 
compare the LGRB rate to the cosmic star formation rate (Sec. l4.lt . 
and the properties of LGRBs hosts to the global properties of galax- 
ies at the same cosmic epoch (Sec. |4.2t . When we study the average 
or median property of the host galaxies we weigh each host by the 
likelihood that it contains a GRB. That allow us to compare with 
observed data since the galaxies which host many GRBs are pro- 
portionally more likely to appear in any given observed sample. We 
also study the typical environment (Sec. 14.3b and evolutionary stage 
(Sec.H3} of LGRB host galaxies. 



4.1 The cosmic star formation rate versus the LGRB rate 

The collapsar model links LGRBs to the evolution of single mas- 
sive stars whose lifetimes are negligible on cosmological scales. If 
no other condition is required for producing a LGRB event, then 
the rate of LGRBs should be an unbiased tracer of the global star 



formation in the Universe (e.g. Totani 1997; Wijers et al. 1998; 


Mao & Mo 


1998; Porciani & Madau 2001; Bromm & Loeb 2002; 


Fynbo ct al 


2006; Price et al. 


2006; Savaelidl200dlTotanietal. 


20061; Prochaska et al. 2007; 


Li 2008b, and references therein). 


Fvnbo et al 


( 2008!) have recently suggested that GRB and Damped 



Lyman- Alpha samples, in contrast with magnitude limited samples, 
provide an almost complete census of z ~ 3 star-forming galaxies. 
We note, however, that the sample used in this study are biased 
against high-metallicity and dusty systems. 



However, both observations and theoretical studies indicate that the 
metallicity of the progenitor star plays an important role in set- 
ting the necessary conditions for a LGRB explosion. In this case, 
the rate of LGRBs is expected to be a biased tracer of the cosmic 
star formation rate. This is demonstrated explicitly in Fig. [2] which 
compares the cosmic star formation rate obtained using all galax- 
ies in the simulation box to that obtained for the three LGRB host 
samples defined in the Sec. [3] 

The sample with no threshold on metallicity (HOST1) traces 
exactly the global star formation rate (solid line in Fig. [2j- This 
is not the case for the two samples with metallicity thresholds 
(HOST2 - dot-dashed line and HOST3 - dashed line). For the 
HOST2 and HOST3 samples, the LGRB rate peaks at higher red- 
shift than the cosmic star formation rate, as a consequence of the 
global decrease of metallicity with increasing redshift. At higher 
redshift, the mean metallicity of the intergalactic medium is lower, 
which implies that a larger fraction of stars form below the metal- 
licity thresholds adopted for HOST2 and HOST3. Therefore the 
deviation of the LGRB rate from the star formation rate decreases 
with increasing redshift. At z ~ 9, the HOST2 sample measures 
about 99 per cent of the global star formation density. The fraction 
decreases to about 30 per cent at z ~ 5 and to only about 10 per 
cent at present. For the HOST3 sample, the bias is even stronger 
because of the lower metallicity threshold: it traces only 30 per 
cent of the global star formation density at z ~ 9, and less than 
10 per cent at z < 5. The results shown in Fig. [2] are in quali- 
tative agreement with recent observational estimates ( Kist ler et al.l 
2008), and with recent theoretical studies also based on the col- 
lapsar model iYoon. Langer & Normanl 120061 ICen & Fan3l2007l. 
iNuza et all 20071 lLapi et alj|2008l) . 

Finally, Fig.f2]also shows that adopting a cosmological model 
compatible with third-year WMAP measurements (right panel), a 
delay is produced in the cosmic star formation rate, due to the delay 



in structure formation. Except for this delay, the predicted trends 
are the same for the WMAP1 and WMAP3 simulations. In the fol- 
lowing, we will only show results obtained using the WMAP3 sim- 
ulation as those obtained for the WMAP1 simulation are very sim- 
ilar. In addition, we will focus only on the two host samples with 
metallicity thresholds (HOST2 and HOST3). 

We remind the reader that we consider as 'hosts' all galaxies 
which can host at least one LGRB event between two simulation 
output. Since galaxies at higher redshift have lower metallicities 
and form stars at higher rates, the rate of LGRBs per host galaxy 
increases rapidly with redshift. The left panel in Fig. [3] shows the 
redshift evolution of the number of LGRBs (solid line) and of host 
galaxies (HOST3 - dashed line). The two vertical lines indicate the 
peaks of the distributions: the number of LGRBs peaks at z ~ 5, 
while the number of host galaxies is maximum at z ~ 2. The right 
panel of Fig. [3] shows the rate of LGRBs per galaxy and per Myr 
computed for the WMAP3 simulation and for the HOST3 sample. 
The predicted rate of LGRBs decreases from ~ 1000 at z ~ 9 
to about 1 Myr -1 galaxy - 1 at z ~ 0, in agreement with calcu- 
lations by iFrver et all dl999l) . Note that the LGRB rate computed 
above is not directly comparable to the observed rate because that 
would require us to take into account many unknown factors like 
the jet angle, an d t o inclu de any possible observational bias (see 
lLapi et al.f2008l : lLHl2008bl) . 

4.2 Physical properties of LGRB host galaxies 

A number of recent papers have studied the physical proper- 
ties of LGRB host galaxies using deep observations covering a 
large wavelength range, both in imaging and in spectroscopy. 
These studies have revealed that LGRB host galaxies are typically 
faint and star forming galaxies, dominated by young and metal- 
poor stellar popul ations (Le Floc'h et alj2003l ; Fruc hter et 
IWainwright etai] |2007|; ISavaglio et ai]|2008l) . In this section, we 
analyse the physical properties of our model host galaxies, and 
compare these with the properties of the average galaxy popula- 
tion, and with observational estimates. 

Fig. [4] shows the K-band rest-frame luminosity function of 
host galaxies for the HOST2 (dot-dashed line) and the HOST3 
(dashed line) samples, compared with the galaxy luminosity func- 
tion measured considering all galaxies in the simulation box, at 
different redshift. At all redshifts, LGRB host galaxies have lumi- 
nosities well below the characteristic luminosity L», in agreement 
with observational measurements. While the total luminosity func- 
tion evolves strongly with the redshift (particularly beyond z ~ 1), 
the number densities and the range of luminosities of LGRB host 
galaxies vary more mildly, due to the fact that at higher redshift a 
larger fraction of the whole galaxy population can host LGRBs. 

Recent observational studies have focused on the stellar mass 
distribution of GRB host galaxies. ICastro Ceron et al" have 
found that the typical stellar mass of host galaxies is smaller than 
the stellar mass of field galaxies at the same redshift. For a sample 
of 30 LGRB hosts, they provide estimates of the stellar mass be- 
tween 10 7 and 10 11 M Q , with a mean value of M» ~ 10 9 7 M . 
About 70 per cent of the host galaxies in their sample have stel- 
lar mass M* < 10 1 01 M@. Similar results have been found by 
ISavaglio et alj d2008). Using a sample of 46 GRB hosts -the largest 
sample so far- they estimate a median stellar mass of 10 9 ' 3 M Q , 
and find that about 83 per cent of the studied systems have stellar 
mass between 10 8 5 and 10 10 3 M . 

In Fig. [5] we compare the galaxy mass distribution for model 
host galaxies from the HOST2 and HOST3 samples to the distribu- 
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Figure 2. Log SFR [Mq yr — 1 Mpc -3 ] as a function of redshift, computed for all galaxies in the simulation box (solid line). Left and right panels correspond 
to the WMAP1 and WMAP3 simulations respectively. The dot-dashed and dashed lines corresponds to the HOST2 and HOST3 samples. The sample with no 
threshold on metallicity (HOST1) traces exactly the global star formation rate measured considering all galaxies in the simulated boxes (solid line). 
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Figure 3. Left panel: redshift distribution of LGRB events (solid line) and of host galaxies (dashed line), for the sample HOST3 (with metallicity threshold 
0.1 Zq). The vertical lines indicate the peaks of the distributions. Right panel: the rate of LGRBs per galaxy per Myr. Both panels correspond to results from 
the WMAP3 simulation. 



tion obtained considering all galaxies in the simulated box. For this 
figure, all galaxies and hosts at all redshifts up to z ~ 9 have been 
used. Fig.[5]shows that LGRB host galaxies have typically low mass 
with a small fraction of them having stellar mass up to ~ 10 11 Mq . 
About 90 per cent of the host galaxies have stellar mass < 10 9 Mq 
and < 10 10 Mq for HOST3 and HOST2 respectively. 

It is well known that the galaxy stellar mass is tightly corre- 
lated with the rest-frame K-band luminosity. ISavaglioe t al. (2008) 
have shown that this relation applies to GRB host galaxies as well, 
but they argue that GRB galaxies have on average higher lumi- 
nosity than "field" galaxies with the same stellar mass, implying a 
lower M*/Lk ratio, as expected for younger galaxies. We compare 
results from our model to observational measurements in Fig. [6] 
The dashed black lin e is the best-fit to the observational data by 
ISavaglio et all d2008f) : logM* = -0.463 x M K - 0.102. The red 
line in Fig. [6] shows the mean luminosity-mass relation obtained 
by using all galaxies in the simulation boxes up to z ~ 9. The 



pink line shows the mean value for host galaxies in the HOST2 
sample, and the blue line corresponds to the mean value obtained 
for the HOST3 sample. To compute the average of Mk we weigh 
each host by the likelihood that it contains a GRB. For this sam- 
ple, we also show th e quartiles of the distribution. We note that 
ISavaglio etal] {2008) adopt a Baldry & Glazebrook IMF for their 
stellar mass estimates, while the model used in this study adopts a 
Chabrier IMF to compute model magnitudes. In order to compare 
model results with observational estimates, we have decreased the 
observed stellar mass by a factor 1.3. Fig.|6]shows that the K-band 
absolute magnitude distribution of simulated GRB host galaxies is 
in good agreement with observations. It also shows that, on aver- 
age, host galaxies have stellar masses which are lower, although 
with a large scatter, than "typical" galaxies with the same mass, in 
agreement with observational findings. 

In Fig. [7] we compare the median colou r of model LGRB host 
galaxies with observational measurements by Savagl io et al. (2008, 
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Figure 4. Luminosity function of the HOST2 (dot-dashed line) and the HOST3 (dashed line) samples, compared to the galaxy luminosity function measured 
using all galaxies in the simulation box (solid line). Different panels are for different redshifts: (top left panel), 1.08 (top right panel), 3.06 (bottom left 
panel), and 6.20 (bottom right panel). In the last two panel the dark bullet shows the characteristic luminosity L*. 
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HOST3 sample (black crosses) compared with the observational 
estimates for the GRB-DLAs studied in Savaglio et al. (2006, 2008 
- blue triangles with error bars a few show the lower and upper- 
branch metallicity solution in Savaglio et al. 2008). In order to com- 
pare with observations we weigh each host with the total number 
of LGRBs. The red line in Fig.[8]shows the median metallicity ob- 
tained using all galaxies in the simulation box. The figure shows 
that the metallicity of model galaxies (both "normal" and hosts) 
does not evolve significantly with redshift. The observational mea- 
surements exhibit a large scatter and have typically large uncertain- 
ties. Within these, model predictions are in relatively good agree- 
ment with observational data. It should be noted that the lack of 
evolution in the gas-phase metallicity of the HOST3 sample is es- 
sentially due to our selection method. In order to enter this sample, 
host galaxies must have young stars with metallicity lower than 
Z ^ O.IZq. In the model used in this study, stars form with the 
metallicity of the gas-phase component, so the adopted selection 
requires host galaxies to have gas-phase metallicity close to the 
adopted threshold (it will be typically higher because the model 
adopts an instantaneous recycling approximation). 



Figure 5. Galaxy mass function for the HOST2 (dot-dashed line) and 
HOST3 (dashed line) samples, compared to the galaxy mass function mea- 
sured using all galaxies in the simulation box (solid line). All galaxies with 
2^9 have been used in this figure. 



shown as black symbols). Model results indicate that GRB galaxies 
are typically bluer than the average galaxy population at the same 
redshift. We note that the observed colours exhibit a quite large 
scatter, probably due to the unknown dust extinction. 

Fig. [8] shows the median gas metallicity evolution for the 



Another important and well studied relation is the luminosity- 
metallicity relation. This has been measured for a very large sample 
of star forming galaxies f rom the Sloan Digital Sky Survey (SDSS) 
by iTremonti et al.1 d2004l) . a nd for smaller samples of galax ies by 
other authors. In particular, iBrown. Kewlev &" Geller (2008) have 
recently suggested a technique to identify extremely metal-poor 
galaxies which share very similar properties (age, m etallicity, star 
forma tion rates) with hosts of LGRBs. The data from lBrown et al.1 
(2008) are plotted in Fig. [9] as blue circ les, tog ether with the ob- 
served relations by Tremonti et alJ (2004b, and bv lRicher & McCalil 
( 1995) for a sample of irregular galaxies. Model results for the 
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Figure 6. Mean K-band absolute magnitude as a function of stellar mass for 
all galaxies in the simulation (red line), for galaxies in the HOST2 sample 
(green), and for galaxies in the HOST3 sample (blue). For the HOST2 and 
HOST3 samples we compute the average of Mr weighing each host with 
the total number of LGRBs. B lack symbols correspond to observational 
data from Savaglio et al. (2008), and the dotted black line is the fit provided 
by the same author (log M* = -0.463 X M K - 0.102). 



Figure 8. Median gas metallicity as a function of redshift for all galaxies 
in the simulation box (red line), and for galaxies from the HOST3 sample 
(dark crosses, error bars show the quartiles of the distributions). For HOST3 
sample we compute the median gas metallicity weighing each host with the 
total number of LGRBs. The blue tria ngles are obs ervational measurements 
for GRB-DLAs h ost galaxies fromlSavagliol j2006h and for the GRB host 
galaxies studied in Sava glio et al.l 120081) . 



CQ 




Figure 7. Median B-K colour as a function of redshift for observed GRB 
hosts (black crosses), all galaxies in the simulation (solid line), galaxies 
from the HOST2 sample (dot-dashed line), and galaxies from the HOST3 
sample (red crosses). For the HOST2 and HOST3 samples we compute the 
median weighing each host with the total number of LGRBs. 



HOST3 sample are plotted as magenta asterisks and lie in the same 
region occupied by the Brown et al. data. Dark crosses show the 
corresponding results for the HOST2 sample. As explained above, 
the adopted selection results in clear metallicity cuts in Fig. [9] Ob- 
servational studies of gas-phase metallicity of GRB hosts could 
therefore provide important information on the metallicity of the 
progenitor stars, although inhomogeneous mixing of metals could 
complicate the interpretation. The few objects in the HOST3 sam- 
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Figure 9. Gas metallicity as function of the B-band luminosity for model 
galaxies with z < 0.1 from the HOST3 (magenta asterisks) and HOST2 
(black crosses) sam ples, compared with observational measurements by 
iBrown et alj j2008l blue circles) for a sample of metal-poor galaxies. The 
blue-dashed and black-dotted li nes show the metallicity-luminosity relation 
measured by Richer & McCall (1995) for a sample of dwarf galaxies, and 
bv lTremonti et alj )2004l) for SPSS star forming galaxies. 
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pies with high metallicity are galaxies with very high star formation 
rate (which coupled with the instantaneous recycling approxima- 
tion, results in quite high metallicities of the inter-stellar medium). 



4.3 The environment of LGRB host galaxies 

The physical environment of LGRB host galaxies can provide 
important information on the origin of LGRBs. The analysis of 
GRB hosts environments is, however, quite difficult from the 
observational viewpoint, given the very low number of identi- 
fied and well studied host galaxies. Only a few observational 
studie s have attempted to a ddress this question in the past few 
years. iGorosabel et alj d2003l) used photometric redshift informa- 
tion in a field of 6' x 6' containing the host galaxy of GRB 
00021 0, and found no obvio us galaxy concentration around the 
host. iBornancini et al.1 d2004h analysed the cross-correlation func- 
tion between host galaxies and surrounding field galaxies using 
VLT and public HST data, and concluded that host galaxies do not 
reside in high density environments. 

The semi-analytic catalogues used in our study provide the 
positions of the LGRB host galaxies, as well as of all other galaxies 
in the simulation box. We have used this information to study the 
auto-correlation function of host galaxies, and the cross-correlation 
function between hosts and all galaxies in the simulation (which we 
will call "normal" galaxies). 

In o rder to compute the tw o-point correlation function, we 
adopt the lLandv & Szalavl J 19931) estimator: 



DD(r) - 2DR(r) 
RR(r) 



(3) 



where DD(r) is the number of galaxy-galaxy pairs at distance r, 
RR(r) is the number of random-random pairs, and DR(r) is the 
number of random-galaxy pairs. 

Fig. [10] shows the auto-correlation function for LGRB host 
galaxies in the HOST2 sample (dashed line) and for normal galax- 
ies (solid line). The number of galaxies in the HOST3 sample 
are too low to compute a reliable correlation function. Results 
are shown for the WMAP3B model, but they are similar for the 
WMAP1 simulation. At low redshift, the auto-correlation function 
of host galaxies is lower than the corresponding function for normal 
galaxies. The difference between the two functions decreases with 
increasing redshift and the two functions almost perfectly overlap 
at z > 4. As explained in Sec. 14.11 this is due to the decrease of 
metallicity at high redshift which implies that an increasing fraction 
of the global galaxy population can host a LGRB event according 
to our selection (see Sec. [3}. 

The two point auto-correlation functions for the HOST2 sam- 
ple and for normal galaxies at z = is repeated in Fig.Qj] together 
with the cross-correlation function between host and normal galax- 
ies (solid line). The cross-correlation function lies in between the 
two auto-correlation functions, suggesting that the probability of 
finding another host near a GRB host is lower than the correspond- 
ing probability of finding a normal galaxy. 

Our results are in qualitative agreement with those found by 
IBornancini et alj d2004h and suggest that LGRB host galaxies tend 
to populate regions with density lower than average. This is not en- 
tirely surprising if one considers that host galaxies are typically 
low-mass star formin g galaxies which pref erentially live in low 
density environments ( Kauffmann et al. 20Q4])- 
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Figure 11. Two-point auto-correlation function for HOST2 galaxies 
(dashed line) and for normal galaxies (dot-dashed line) in the WMAP3B 
model. The solid line shows the cross-correlation function between host 
and normal galaxies. 



4.4 Descendants of high-z LGRB host galaxies 

From the theoretical point of view, it is interesting to ask which are 
the 'descendants' of high-redshift LGRB host galaxies. Do they 
preferentially end up in massive haloes? What are the typical mor- 
phology, colour, mass and metallicity of the descendants? While 
this is a very difficult (if not impossible) question to address ob- 
servationally, it can be easily addressed with the available semi- 
analytic catalogues, which contain the full merger tree information 
for all galaxies in the simulation box. 

In this section, we use this information to study the fate of 
LBRG host galaxies selected at z ~ 4 (hereafter progenitors) in 
the observed mass-metallicity and colour-magnitude planes, and 
the distribution of host halo virial mass of the descendant galax- 
ies. For simplicity, we only show results for our HOST3 sample 
which contains a lower number of host galaxies than the HOST2 
sample. Results are, however, similar for HOST2. 

The top-left panel in Fig. [J_2] shows the mass-metallicity rela- 
tion for our HOST3 sample at z = 4.18. The other panels show 
the location in the same plane of the descendant galaxies, down to 
z = 0. Galaxies are colour-coded as functions of their morpholog- 
ical types which are assigned by the ratio of the total bulge mass 
to the total stellar mass, B/T. Objects with B/T > 0.5 are classi- 
fied as 'ellipticals' and are shown as red triangles, black diamonds 
represent 'spirals' (B/T < 0.5). 

Fig. [J_2] shows that the LGRB host galaxies at redshift ~ 4 
are low-mass and low-metallicity galaxies (as discussed in the pre- 
vious sections), and that the great majority of them do not have a 
significant bulge component. At later time, galaxies grow in mass 
and their gas-phase metallicity increases. Fig|12lshows that the in- 
crease of the gas-phase metallicity is very rapid: a relatively large 
fraction of descenda nts at z ~ 3 already have solar metallicity 
[12+log (O/H) ~8.7, lAllende Prieto et alj|200ill and most of the 
descendants at z ~ 2 have super-solar metallicity. This efficient 
enrichment of the interstellar medium is due to the fact that the 
semi-analytic model used in our study assumes a perfect mixing ef- 
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Figure 10. Two-point auto-correlation function for LGRB host galaxies from the HOST2 sample (dashed line) and for normal galaxies (solid line), at different 
redshifts, for WMAP3B cosmology. 



ficiency of the metals formed by new stars (De Lucia et ai] |2004 . 
The stellar mass of the descendants evolves more slowly, with very 
few objects jumping to masses larger than 10 11 Mq, probably as 
a consequence of major mergers. At z = 0, all descendants of 
LGRBs at z ~ 4 have relatively high gas-phase metallicities and 
cover all the mass range between a few times 10 9 and 10 12 Mq. 
Interestingly, a large fraction of these (about 66 per cent) have a 
dominant bulge component. 

Fig. Q~3] shows the location in the colour-magnitude diagram 
of LGRB host galaxies at z = 4.18 (top left panel) and of all their 
descendants at later time. LGRB hosts selected at z ~ 4 have very 
blue colours (R-K < 1) and relatively faint magnitudes. The de- 
scendants of these galaxies become progressively redder (they all 
have R-K > 2 at z — 0) and cover a wider range of magnitudes 
at low redshift. A few objects become very luminous and very red 
already at z ~ 3 (these are the same objects which appear very 
massive and metal rich at the same redshift in Fig.l 121. 

Finally, in Fig.[l4]we show the parent halo mass of all descen- 
dant of LGRB host galaxies at z ~ 4 versus the halo mass of the 
host. The blue dashed line in each panel shows the one-to-one rela- 
tion and is plotted to guide eyes. As expected, descendant galaxies 
reside in larger and larger haloes as time goes on. At z ~ 2, the 
majority of the descendant galaxies still reside in haloes of mass 
~ 10 11 M e but a few of them are in relatively massive haloes 
(~ 10 14 Mq) and are satellite galaxies of a very massive cluster 
at z = 0. In the local Universe, the descendants of high-z LGRB 
hosts reside in haloes of different mass but most of them still reside 



in haloes with mass between 10 and 10 Mq, in agreement with 
the results in Sec. 14. 31 



5 DISCUSSION AND CONCLUSIONS 

In this work, we have studied the physical and environmental prop- 
erties of galaxies hosting Long Gamma-Ray Bursts (LGRBs) in the 
context of a hierarchical model of galaxy formation. In order to se- 
lect host galaxies, we have adopted the collapsar model and used 
information available from a semi-analytic model of galaxy forma- 
tion coupled to high-resolution cosmological simulations (Wang et 
al. 2008). 

By imposing different metallicity constraints on the progeni- 
tor stars, we have constructed three host galaxies samples: HOST1 
is built without any cut on the metallicity of progenitor stars of 
GRBSs, while the HOST2 and HOST3 samples are constructed by 
selecting galaxies with progenitor stars of metallicity lower than 
0.3Zq and O.IZq respectively. 

A number of recent studies have adopted a sim ilar but not 
identic al approach to study the host galaxies of LGRBs. iNuza et al.l 
J2007I) developed a Monte Carlo code to identify hosts of LGRBs 
within a cosmological hydrodynamical simulations. Their analysis 
was also based on the collapsar model and limited to 2 ^ 3. The 
simulation used in Nuza et al. was relatively small (10/i _1 Mpc) 
and therefore did not allow a detailed investigation of the envi- 
ronmental properties of LGRB hosts. As confirmed by our study, 
they pointed out that if LGRBs are generated by the death of mas- 



© 2009 RAS, MNRAS 000,[TV?? 



10 Campisi M.A. etal. 




Figure 12. Mass-metallicity relation for HOST3 galaxies at redshift 4.18 (top left panel), and for their descendants at lower redshifts in the WMAP3B model. 
The colour coding indicates galaxies with different morphologies (re d trian gles for ellipticals and black diamonds for spirals). The solid black line in the 
bottom right panel shows the best fit relation found by Tremonti et al. ( 2004). The horizontal dashed line corresponds to the adopted solar abundance. 





Figure 13. R — K colour as a function of the K-band magnitude for HOST3 galaxies galaxies at z ~ 4 (top-left panel), and for their descendants at lower 
redshifts in the WMAP3B model. Red triangles are for elliptical galaxies while black diamonds are for spirals. 
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Figure 14. Parent halo mass of all descendants of HOST3 galaxies selected at z = 4. 18, as a function of the halo mass of the host galaxies in the WMAP3B 
case. Red symbols are used for satellite galaxies and black symbols for central galaxies. 



sive young stars of metallicity lower than a certain value, they do 
not provide a good tracer of the cosmic star formation history (see 
Sec.|4jT]and later). 

lLapi et all d2008l) used the cosmological star formation rate 
below a critical metallicity to estimate the event rate of LGRBs. To 
this purpose, they employed the galaxy formation model presented 
in Granato et al. (2004). Lapi et al. find that their predicted number 
counts of LGRBs agrees well with the bright SWIFT data, without 
the need for an intrinsic luminosity evolution. They find that host 
galaxies are dominated by young stellar populations, are gas rich 
and metal-poor. The model adopted in Lapi et al. (2008) does not 
follow galactic disks nor does it consider mergers between galaxies 
or between haloes. In addition, the model does not provide spatial 
information of model galaxies. 

In another recent study, ICourtv et all (2007) used N- 
body/Eulerian hydrodynamic simulations, and identified GRB host 
galaxies with those having SFR and SFR-to-luminosity similar to 
those of 10 observed GRB hosts in the redshift range 0.43 < z < 
2.03. They found that the host galaxies have low stellar masses and 
low mass-to-light ratios, are young and bluer than typical galax- 
ies at the same cosmic epochs. Their identification of simulated 
galaxies with observed hosts is limited by uncertainties on the ob- 
servational estimates of the SFR and luminosity of host galaxies. In 
addition, their simulated volume is relatively small (32 Mpc) 
and the analysis is limited to relatively low redshifts. 

Compared to previous works, our study uses a much larger 
simulated volume (a box of 125 hT 1 Mpc on a side - see Sec.[2j 
that allows us to study the environmental properties of model host 
galaxies. In addition, the use of two different cosmological models 
allows us to analyse the dependence of results on cosmology. The 
semi-analytic model employed in our work has been studied in a 
number of previous papers, and it has been shown to successfully 



reproduce a number of observational results for the global galaxy 
population in the local Universe and at higher redshift (see Sec.[2]l. 

In agreement with previous work, and as expected due to the 
global increase of the ISM metallicity with decreasing redshift, we 
find that when assuming a metallicity threshold for progenitor stars 
of LGRBs, they do not represent a perfect tracer of the cosmic star 
formation history. The bias is stronger as the metallicity threshold 
assumed is lowered. At higher redshift, the cold-phase metallicity 
of model galaxies is lower, and galaxies form stars at higher rates. 
As a consequence, there is a higher rate of LGRBs per galaxy and 
the host galaxy population includes a larger fraction of the global 
galaxy population, with the two populations sharing very similar 
physical properties. 

At lower redshift, the host galaxy population is dominated 
by galaxies with low-masses, relatively young ages, blue colours, 
and luminosity below L», in qualitative agreement with observa- 
tional measurements. We note, however, that while ~ 90 per cent 
of the galaxies in our simulations have stellar masses lower than 
~ 10 10 Mq, only ~ 0.3 per cent of them are in the HOST3 sample 
and ~ 13 per cent in the HOST2 sample. 

The metal content and metallicity evolution of host galaxies 
also appear to be in agreement with observational estimates. Since 
in the model adopted in this study stars form with the metallicity 
of the cold-phase component at the time of the star formation, host 
galaxies have average metallicity lower or close to the metallic- 
ity threshold adopted. Metallicity studies of LGRB host galaxies 
could therefore provide important constraints on the adopted col- 
lapsar model. It is important to note, however, that inhomogeneous 
mixing of metals (which is not considered in the model adopted 
here) could significantly complicate this picture. The metal content 
of host galaxies does not significantly evolve with redshift, while 
the overall galaxy population exhibit some weak evolution. 
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Taking advantage of the spatial information provided by the 
semi-analytic model used in our work, we study the clustering 
properties of model host galaxies and compare them to the corre- 
sponding properties of the global galaxy population. Since a larger 
fraction of the whole galaxy population can host LGRBs at high 
redshift, the clustering properties of host galaxies do not differ sig- 
nificantly from the clustering properties of 'normal' galaxies. At 
lower redshift, the host galaxies are significantly less clustered than 
normal galaxies, as expected due to their physical properties. Inter- 
estingly, we find that the cross-correlation function between host 
and normal galaxies lies in between the the auto-correlation func- 
tions of host galaxies and normal galaxies, suggesting that the prob- 
ability of finding another host galaxy nearby a GRB host is lower 
than the corresponding probability of finding a normal galaxy. This 
implies that LGRBs reside in regions with density lower than av- 
erage. When larger samples of LGRB host galaxies will become 
available, it will be possible to test our clustering predictions. 

Using the available semi-analytic catalogues, which contain 
the full merger tree information for all galaxies in the simulation 
box, we find that LGRB host galaxies at redshift ~ 4 (which have 
low stellar mass, low metallicity and no significant bulge compo- 
nent) evolve into relatively massive, red galaxies at redshift zero. 
While the descendants of high-z LGRBs reside in haloes of differ- 
ent mass in the local Universe, most of them still sit in haloes with 
mass between 10 12 and 10 13 M . 
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